[1] "C:/Users/happy/Dropbox/Penn State/Dissertation/splits/Joo_Mukherjee_Replication"
RStudioGD 
        2 
RStudioGD 
        2 
RStudioGD 
        2 
RStudioGD 
        2 
RStudioGD 
        2 
RStudioGD 
        2 
RStudioGD 
        2 

	Welch Two Sample t-test

data:  dat$max_year by dat$cent
t = -0.0052847, df = 117.78, p-value = 0.9958
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -5.395070  5.366352
sample estimates:
mean in group 1 mean in group 3 
       18.29333        18.30769 


	Welch Two Sample t-test

data:  dat$max_year by dat$cent
t = -1.3634, df = 159.85, p-value = 0.1747
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -7.970238  1.460130
sample estimates:
mean in group 1 mean in group 2 
       18.29333        21.54839 


	Welch Two Sample t-test

data:  dat$max_year by dat$cent
t = 1.3079, df = 110.77, p-value = 0.1936
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -1.669151  8.150540
sample estimates:
mean in group 2 mean in group 3 
       21.54839        18.30769 


	3-sample test for equality of proportions without continuity correction

data:  prop.hi out of trial
X-squared = 0.74311, df = 2, p-value = 0.6897
alternative hypothesis: two.sided
sample estimates:
   prop 1    prop 2    prop 3 
0.1538462 0.2197802 0.1860465 


	3-sample test for equality of proportions without continuity correction

data:  prop.low out of trial
X-squared = 2.3459, df = 2, p-value = 0.3095
alternative hypothesis: two.sided
sample estimates:
   prop 1    prop 2    prop 3 
0.4230769 0.2912088 0.2558140 


	3-sample test for equality of proportions without continuity correction

data:  prop.mid out of trial
X-squared = 1.249, df = 2, p-value = 0.5355
alternative hypothesis: two.sided
sample estimates:
   prop 1    prop 2    prop 3 
0.4230769 0.4890110 0.5581395 

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi, data = jm)

  n= 3725, number of events= 144 
   (25 observations deleted due to missingness)

          coef exp(coef) se(coef)     z Pr(>|z|)   
cent_bi 0.5441    1.7231   0.1849 2.943  0.00325 **
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

        exp(coef) exp(-coef) lower .95 upper .95
cent_bi     1.723     0.5804     1.199     2.475

Concordance= 0.525  (se = 0.019 )
Likelihood ratio test= 7.96  on 1 df,   p=0.005
Wald test            = 8.66  on 1 df,   p=0.003
Score (logrank) test = 8.88  on 1 df,   p=0.003

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr, data = jm)

  n= 3488, number of events= 140 
   (262 observations deleted due to missingness)

            coef exp(coef) se(coef)      z Pr(>|z|)   
cent_bi   0.5493    1.7321   0.1904  2.885  0.00391 **
mobilize  0.1836    1.2016   0.1425  1.289  0.19750   
terr     -0.3643    0.6946   0.1862 -1.957  0.05037 . 
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

         exp(coef) exp(-coef) lower .95 upper .95
cent_bi     1.7321     0.5773    1.1927     2.516
mobilize    1.2016     0.8322    0.9088     1.589
terr        0.6946     1.4396    0.4823     1.001

Concordance= 0.54  (se = 0.025 )
Likelihood ratio test= 11.58  on 3 df,   p=0.009
Wald test            = 11.99  on 3 df,   p=0.007
Score (logrank) test = 12.18  on 3 df,   p=0.007

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire, data = jm)

  n= 3481, number of events= 140 
   (269 observations deleted due to missingness)

               coef exp(coef) se(coef)      z Pr(>|z|)   
cent_bi      0.4885    1.6298   0.1911  2.557  0.01057 * 
mobilize     0.1070    1.1129   0.1472  0.727  0.46743   
terr        -0.4051    0.6669   0.1865 -2.172  0.02985 * 
negotiation  0.5165    1.6761   0.1940  2.663  0.00775 **
ceasefire    1.0905    2.9759   0.4216  2.587  0.00969 **
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

            exp(coef) exp(-coef) lower .95 upper .95
cent_bi        1.6298     0.6136    1.1207    2.3702
mobilize       1.1129     0.8986    0.8340    1.4851
terr           0.6669     1.4995    0.4627    0.9612
negotiation    1.6761     0.5966    1.1461    2.4513
ceasefire      2.9759     0.3360    1.3024    6.7997

Concordance= 0.605  (se = 0.026 )
Likelihood ratio test= 24.15  on 5 df,   p=2e-04
Wald test            = 28.23  on 5 df,   p=3e-05
Score (logrank) test = 29.85  on 5 df,   p=2e-05

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2, 
    data = jm)

  n= 2590, number of events= 120 
   (1160 observations deleted due to missingness)

                   coef  exp(coef)   se(coef)      z Pr(>|z|)    
cent_bi       0.4920553  1.6356745  0.2120456  2.321 0.020313 *  
mobilize      0.2950072  1.3431360  0.1610421  1.832 0.066972 .  
terr         -0.6963386  0.4984068  0.2054148 -3.390 0.000699 ***
negotiation   0.3769187  1.4577858  0.2038009  1.849 0.064394 .  
ceasefire     0.6974000  2.0085237  0.4402755  1.584 0.113192    
independence -1.0088887  0.3646240  0.2300596 -4.385 1.16e-05 ***
logbd         0.0323538  1.0328829  0.0384169  0.842 0.399689    
logbd2       -0.0007152  0.9992851  0.0044482 -0.161 0.872272    
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

             exp(coef) exp(-coef) lower .95 upper .95
cent_bi         1.6357     0.6114    1.0795    2.4785
mobilize        1.3431     0.7445    0.9796    1.8416
terr            0.4984     2.0064    0.3332    0.7455
negotiation     1.4578     0.6860    0.9777    2.1735
ceasefire       2.0085     0.4979    0.8474    4.7604
independence    0.3646     2.7426    0.2323    0.5724
logbd           1.0329     0.9682    0.9580    1.1137
logbd2          0.9993     1.0007    0.9906    1.0080

Concordance= 0.714  (se = 0.024 )
Likelihood ratio test= 66.36  on 8 df,   p=3e-11
Wald test            = 68.12  on 8 df,   p=1e-11
Score (logrank) test = 74.14  on 8 df,   p=7e-13

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2 + 
    state_svac_reb_lag + log_os_reb + log_os_gov, data = jm)

  n= 2391, number of events= 118 
   (1359 observations deleted due to missingness)

                        coef exp(coef)  se(coef)      z Pr(>|z|)    
cent_bi             0.467661  1.596256  0.217147  2.154  0.03127 *  
mobilize            0.300123  1.350024  0.165055  1.818  0.06902 .  
terr               -0.708063  0.492597  0.215699 -3.283  0.00103 ** 
negotiation         0.376467  1.457127  0.207012  1.819  0.06898 .  
ceasefire           0.791768  2.207296  0.446600  1.773  0.07625 .  
independence       -0.988190  0.372250  0.230886 -4.280 1.87e-05 ***
logbd               0.002957  1.002961  0.044400  0.067  0.94691    
logbd2             -0.003433  0.996573  0.004893 -0.702  0.48285    
state_svac_reb_lag  0.124468  1.132546  0.187781  0.663  0.50744    
log_os_reb          0.021000  1.021222  0.015956  1.316  0.18814    
log_os_gov          0.009753  1.009801  0.011318  0.862  0.38881    
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

                   exp(coef) exp(-coef) lower .95 upper .95
cent_bi               1.5963     0.6265    1.0430    2.4431
mobilize              1.3500     0.7407    0.9769    1.8657
terr                  0.4926     2.0301    0.3228    0.7518
negotiation           1.4571     0.6863    0.9712    2.1863
ceasefire             2.2073     0.4530    0.9198    5.2967
independence          0.3722     2.6864    0.2368    0.5853
logbd                 1.0030     0.9970    0.9194    1.0942
logbd2                0.9966     1.0034    0.9871    1.0062
state_svac_reb_lag    1.1325     0.8830    0.7838    1.6364
log_os_reb            1.0212     0.9792    0.9898    1.0537
log_os_gov            1.0098     0.9903    0.9876    1.0325

Concordance= 0.727  (se = 0.024 )
Likelihood ratio test= 70.82  on 11 df,   p=9e-11
Wald test            = 74.63  on 11 df,   p=2e-11
Score (logrank) test = 82.82  on 11 df,   p=4e-13

Call:
coxph(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2 + 
    state_svac_reb_lag + log_os_reb + log_os_gov + ngroup_country + 
    polwing, data = jm)

  n= 2351, number of events= 116 
   (1399 observations deleted due to missingness)

                        coef exp(coef)  se(coef)      z Pr(>|z|)    
cent_bi             0.398296  1.489285  0.220088  1.810 0.070340 .  
mobilize            0.161776  1.175597  0.176026  0.919 0.358073    
terr               -0.620156  0.537861  0.216227 -2.868 0.004130 ** 
negotiation         0.346567  1.414204  0.211200  1.641 0.100809    
ceasefire           0.833694  2.301807  0.447874  1.861 0.062681 .  
independence       -0.829491  0.436271  0.235527 -3.522 0.000429 ***
logbd              -0.012390  0.987686  0.046327 -0.267 0.789119    
logbd2             -0.005172  0.994842  0.005075 -1.019 0.308196    
state_svac_reb_lag  0.068862  1.071288  0.188962  0.364 0.715543    
log_os_reb          0.014991  1.015104  0.016598  0.903 0.366409    
log_os_gov          0.021272  1.021500  0.011807  1.802 0.071604 .  
ngroup_country     -0.077580  0.925353  0.028090 -2.762 0.005747 ** 
polwing            -0.220606  0.802032  0.208828 -1.056 0.290784    
---
Signif. codes:  0 �***� 0.001 �**� 0.01 �*� 0.05 �.� 0.1 � � 1

                   exp(coef) exp(-coef) lower .95 upper .95
cent_bi               1.4893     0.6715    0.9675    2.2925
mobilize              1.1756     0.8506    0.8326    1.6599
terr                  0.5379     1.8592    0.3521    0.8217
negotiation           1.4142     0.7071    0.9348    2.1394
ceasefire             2.3018     0.4344    0.9568    5.5373
independence          0.4363     2.2922    0.2750    0.6922
logbd                 0.9877     1.0125    0.9020    1.0816
logbd2                0.9948     1.0052    0.9850    1.0048
state_svac_reb_lag    1.0713     0.9335    0.7397    1.5515
log_os_reb            1.0151     0.9851    0.9826    1.0487
log_os_gov            1.0215     0.9790    0.9981    1.0454
ngroup_country        0.9254     1.0807    0.8758    0.9777
polwing               0.8020     1.2468    0.5326    1.2077

Concordance= 0.754  (se = 0.023 )
Likelihood ratio test= 78.66  on 13 df,   p=2e-11
Wald test            = 81.13  on 13 df,   p=7e-12
Score (logrank) test = 89.81  on 13 df,   p=2e-13


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi, data = jm, 
    dist = "weibull")
              Value Std. Error     z      p
(Intercept)  5.9468     0.2006 29.64 <2e-16
cent_bi     -0.5016     0.1663 -3.02 0.0026
Log(scale)  -0.1184     0.0676 -1.75 0.0798

Scale= 0.888 

Weibull distribution
Loglik(model)= -1025.1   Loglik(intercept only)= -1029.4
	Chisq= 8.54 on 1 degrees of freedom, p= 0.0035 
Number of Newton-Raphson Iterations: 9 
n=3725 (25 observations deleted due to missingness)


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr, data = jm, dist = "weibull")
              Value Std. Error     z      p
(Intercept)  5.9980     0.2842 21.11 <2e-16
cent_bi     -0.5029     0.1688 -2.98 0.0029
mobilize    -0.1470     0.1247 -1.18 0.2386
terr         0.3280     0.1626  2.02 0.0436
Log(scale)  -0.1287     0.0688 -1.87 0.0613

Scale= 0.879 

Weibull distribution
Loglik(model)= -991.4   Loglik(intercept only)= -997.5
	Chisq= 12.23 on 3 degrees of freedom, p= 0.0066 
Number of Newton-Raphson Iterations: 8 
n=3488 (262 observations deleted due to missingness)


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire, data = jm, dist = "weibull")
              Value Std. Error     z      p
(Intercept)  6.0104     0.2876 20.90 <2e-16
cent_bi     -0.4474     0.1699 -2.63 0.0085
mobilize    -0.0791     0.1290 -0.61 0.5398
terr         0.3633     0.1634  2.22 0.0262
negotiation -0.4470     0.1745 -2.56 0.0104
ceasefire   -0.9622     0.3765 -2.56 0.0106
Log(scale)  -0.1258     0.0689 -1.82 0.0680

Scale= 0.882 

Weibull distribution
Loglik(model)= -985.2   Loglik(intercept only)= -997.5
	Chisq= 24.51 on 5 degrees of freedom, p= 0.00017 
Number of Newton-Raphson Iterations: 8 
n=3481 (269 observations deleted due to missingness)


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2, 
    data = jm, dist = "weibull")
                 Value Std. Error     z       p
(Intercept)   5.526490   0.387292 14.27 < 2e-16
cent_bi      -0.436600   0.186472 -2.34 0.01921
mobilize     -0.232527   0.139843 -1.66 0.09636
terr          0.609211   0.175656  3.47 0.00052
negotiation  -0.318232   0.178523 -1.78 0.07465
ceasefire    -0.616403   0.383770 -1.61 0.10823
independence  0.890307   0.208273  4.27 1.9e-05
logbd        -0.030226   0.033031 -0.92 0.36015
logbd2        0.000369   0.003843  0.10 0.92341
Log(scale)   -0.145060   0.074092 -1.96 0.05025

Scale= 0.865 

Weibull distribution
Loglik(model)= -813.9   Loglik(intercept only)= -847.9
	Chisq= 67.89 on 8 degrees of freedom, p= 1.3e-11 
Number of Newton-Raphson Iterations: 8 
n=2590 (1160 observations deleted due to missingness)


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2 + 
    state_svac_reb_lag + log_os_reb + log_os_gov, data = jm, 
    dist = "weibull")
                      Value Std. Error     z       p
(Intercept)         5.11445    0.46962 10.89 < 2e-16
cent_bi            -0.41253    0.19047 -2.17 0.03033
mobilize           -0.23389    0.14334 -1.63 0.10274
terr                0.62231    0.18622  3.34 0.00083
negotiation        -0.31432    0.18157 -1.73 0.08343
ceasefire          -0.69596    0.38924 -1.79 0.07378
independence        0.86144    0.20902  4.12 3.8e-05
logbd              -0.00664    0.03817 -0.17 0.86192
logbd2              0.00256    0.00422  0.61 0.54338
state_svac_reb_lag -0.10003    0.16146 -0.62 0.53556
log_os_reb         -0.01828    0.01376 -1.33 0.18391
log_os_gov         -0.00828    0.00976 -0.85 0.39634
Log(scale)         -0.14574    0.07521 -1.94 0.05266

Scale= 0.864 

Weibull distribution
Loglik(model)= -795.4   Loglik(intercept only)= -831.5
	Chisq= 72.24 on 11 degrees of freedom, p= 4.6e-11 
Number of Newton-Raphson Iterations: 7 
n=2391 (1359 observations deleted due to missingness)


Call:
survreg(formula = Surv(dur_exist, split) ~ cent_bi + mobilize + 
    terr + negotiation + ceasefire + independence + logbd + logbd2 + 
    state_svac_reb_lag + log_os_reb + log_os_gov + ngroup_country + 
    polwing, data = jm, dist = "weibull")
                      Value Std. Error     z       p
(Intercept)         4.38899    0.50077  8.76 < 2e-16
cent_bi            -0.34739    0.18980 -1.83 0.06720
mobilize           -0.10323    0.15093 -0.68 0.49399
terr                0.53074    0.18534  2.86 0.00419
negotiation        -0.28990    0.18245 -1.59 0.11207
ceasefire          -0.71790    0.38530 -1.86 0.06243
independence        0.71640    0.20690  3.46 0.00054
logbd               0.00573    0.03924  0.15 0.88387
logbd2              0.00394    0.00432  0.91 0.36223
state_svac_reb_lag -0.05103    0.16001 -0.32 0.74978
log_os_reb         -0.01216    0.01416 -0.86 0.39063
log_os_gov         -0.01836    0.01008 -1.82 0.06866
ngroup_country      0.06627    0.02433  2.72 0.00645
polwing             0.19506    0.17559  1.11 0.26661
Log(scale)         -0.15970    0.07596 -2.10 0.03552

Scale= 0.852 

Weibull distribution
Loglik(model)= -778.5   Loglik(intercept only)= -818.5
	Chisq= 79.95 on 13 degrees of freedom, p= 1.1e-11 
Number of Newton-Raphson Iterations: 7 
n=2351 (1399 observations deleted due to missingness)

RStudioGD 
        2 
